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ABSTRACT 

This  report  gives  the  results  of  calculations  of  gamma-ray  dose  rates  in  two-legged 
rectangular  concrete  ducts.  The  calculations  were  performed  on  an  IBM  1620,  They  are 
based  on  the  differential  dose  albedo  and  include  multiple  scattering  effects.  The  results 
of  the  calculations  are  compared  with  experimental  data  for  ducts  whose  widths  vary  from 
11  inches  to  6  feet,  using  Co-60,  Cs-137,  Au-198,  and  Na-24  gamma-ray  sources.  The 
calculated  dose  rates  agree  to  within  ±30  percent  for  all  ducts  and  all  sources,  except 
for  small  ducts  with  Cs-137  sources,  ducts  with  very  short  first  legs  (L^/W  s  1.33),  and 
for  Au-198  sources. 


Qualified  ra(|ua*tar«  may  obtain  eoploi  of  thit  ropert  from  DOC. 
Tho  Loborotory  invito*  comment  on  thi*  report,  particularly  on  the 
reiulti  obtained  by  thoie  who  hove  opplleo  the  Information. 


INTRODUCTION 


The  ability  to  calculate  dose  rates  in  two-legged  ducts  would  be  of  great 
assistance  to  engineers  designing  shelter  entranceways.  If  the  calculations  were 
based  on  basic  scattering  principles  and  were  verified  experimentally,  it  would 
give  great  insight  into  the  phenomena  taking  place  in  radiation  attenuation  In  a 
duct. 


The  LeDoux-Chilton  method^  was  developed  on  the  albedo  concept,  where  the 
dose  contributions  from  the  different  scattering  areas  in  o  duct  were  calculated  using 
the  differential  dose  albedo.  However,  at  the  time  of  the  development  of  this  com¬ 
putational  technique,  values  for  differential  dose  albedo  were  not  known.  Values 
for  differential  dose  albedos  for  various  gamma-ray  energies  and  entrance  and  exit 
angles  have  since  been  calculated  by  Technical  Operations,  using  the  Monte  Carlo 
method. 2  Using  the  Technical  Operations  data,  Chilton  and  Huddleston 3  developed 
o  semiempirical  equation  from  which  the  differential  dose  albedo  can  be  calculated 
for  any  energy  and  any  entrance  and  exit  angle.  Ingold^  used  this  equation  to  cal¬ 
culate  build-up  in  the  first  leg  and  obtained  good  agreement  with  experimental  results. 

With  an  equation  for  albedo,  and  the  equations  for  corner  lip  inscattering  and 
penetration  effects,  It  became  feasible  to  calculate  theoretical  dose  rates  in  concrete 
ducts  with  an  electronic  computer.  For  these  calculations,  four  programs  were  written 
for  the  IBM  1620:  (1)  a  program  which  performs  the  calculations  devised  In  Reference  1, 
(2)  a  program  which  uses  small  incremental  areas  for  calculating  the  dose  rates  from  the 
primary  scattering  areas,  (3)  a  program  which  calculates  the  dose  rates  from  multiple 
scattering  in  the  first  and  second  legs  of  the  duct,  and  (4)  a  program  which  calculates 
the  dose  rates  due  to  multiple  scattering  of  gamma  rays  that  are  deflected  (inscattered) 
by  the  corner  lip.  These  programs  therefore  Include  all  processes  which  involve  one 
or  two  Interactions;  that  is;  one  backscatter,  one  inscatter,  one  backscatter  and  one 
penetration,  two  bockscatters,  and  one  backscatter  and  one  inscatter. 

In  the  study  that  follows,  the  method  developed  for  calculating  the  differential 
dose  albedo  is  presented,  then  the  four  computer  programs  are  outlined,  and  the  results 
of  the  calculations  ore  compared  with  experimental  measurements.  The  study  was 
sponsored  by  the  Defense  Atomic  Support  Agency  through  the  Bureau  of  Yards  and 
Docks. 


CALCULATION  OF  DIFFERENTIAL  DOSE  ALBEDO 


The  dose  rate  at  the  detector  from  the  scattering  area,  A  (Figure  1),  is  given 
by 


A  COS  6a 

- n - -  0) 

^1  ■’2 


where  a(EQ,  0^,  0,  <p)  =  the  differential  dose  albedo 

A  =  the  area  of  the  scattering  surface 
Dq  =  the  dose  rate  at  one  unit  length  from  the  source 
Eq  =  the  initial  energy  of  the  gamma  rays  from  the  source 


Values  for  o(Eo,  0q,  0,  (p)  were  calculated  for  various  energies  and  entrance 
and  exit  angles  by  Technical  Operations,  Inc.,  using  the  Monte  Carlo  method. 2 


Figure  1,  Scattering  of  gamma  rays  from  surface. 
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The  equation  developed  by  Chilton  and  Huddleston^  to  express  the  albedo  for 
a  given  energy  is 


a(EQ,  Sq,  6,  (p) 


26 


C(EQ)K(e^)10"'-  +  C(EQ) 


1  + 


cos  6q 
cos  9 


(2) 


where  C(Eo)  and  C'(Eq)  are  constants  for  a  given  energy,  <(65)  is  the  Klein-Nishina 
differentlol  energy  scattering  coefficient,  and  65  is  the  angle  through  which  the 
radiation  Is  scattered  and  is  given  by 


cos0^  =  sin  6q  sin  6  cos(p  -  cos  Bq  cos  6 


Chilton  and  Huddleston  found  C(Eq)  and  C'(Eq)  by  the  least-squares  fit  of  the 
Technical  Operations  data  by  both  a  nonweighted  fit  and  a  weighted  fit.  The  values 
of  C  and  C  found  by  the  weighted  fit  are  plotted  in  Figure  2.  For  calculation  on  the 
computer,  the  curves  of  Figure  2  were  fitted  by  the  equations: 

C  =  expC-2.921  +  0.6805  InE  +  0.0111  (InE)^  -  0.04131  (In E)®]  (3) 

C  =  exp  [-5. 89  +  0.275  (E  -3.25)^]  0.  Imev  s  E  £  Imev 

C  =  exp  [-4. 86  +  0. 36(E  -  2)^]  Imev  &  E  s  2mev  (4) 

C  =  exp  [-4. 83  -  0. 013E]  2mev  <  E  s;  10  mev 

PROGRAM  1.  LEDOUX-CHILTON  CALCULATIONS 

For  the  calculations  based  on  the  LeDoux-Chilton  method,  ^  the  duct  is  divided 
into  the  areas  shown  in  Figure  3.  A] ,  A2,  and  A3  +  A4  are  referred  to  as  the  primary 
scattering  areas,  those  areas  In  view  of  both  the  source  and  the  detector.  A5 ,  A^ , 

Ay,  and  Ag  are  those  areas  for  which  radiation  must  penetrate  the  corner  lip  before 
reaching  the  detector. 
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Figure  3.  Duct  geometi^  showing  scattering  areas. 
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The  dose  rate  from  any  of  the  scattering  areas  is  calculated  from  Equation  1, 
where  r^ ,  r2 ,  6q,  9,  and  cp,  shown  in  Figure  1,  are  measured  at  the  center  of  the 
scattering  area.  The  albedo,  a(EQ,  9q,  0,  <p),  is  calculated  from  Equation  2. 

C(Eq)  and  C'(Eq)  in  Equation  2  are  calculated  from  Equations  3  and  4. 

For  areas  A^,  A^,  Ay,  and  Ag  the  size  of  the  scattering  area  is  found  by 
assuming  that  all  gamma  rays  traveling  a  distance  of  one  mean  free  path  or  less 
through  the  corner  lip  reach  the  detector,  and  none  traveling  more  than  one  mean 
free  path  reach  the  detector.  The  validity  of  this  assumption  is  shown  in  Reference  1. 

The  mean  free  paths,  or  RL's  (relaxation  lengths),  used  are  the  reciprocals  of 
the  energy-absorption  coefficients,  as  It  is  assumed  that  gamma  rays  scattered  but 
not  absorbed  are  scattered  through  a  small  angle.  These  mean  free  paths,  found  by 
taking  the  reciprocal  of  the  energy-absorption  coefficient  for  air  given  by  Rockwell,^ 
are  plotted  in  Figure  4.  For  calculation  on  the  computer,  the  curve  of  Figure  4  was 
fitted  by  the  equations: 


RL  =  4.65E 

RL  =  5.98 

RL  =  5.98E 


-0.24 


0.269 


E  <  0. 3mev 

O.Smev  £  E  s  l.Omev 

E  >  l.Omev 


0.1  0.2  0.3  0.4  0.5  1  2 

Gaimno-Roy  Energy  (mov) 

Figure  4.  Energy-absorption  mean  free  path  for  145~lb/ft^  air-equivalent 
concrete. 
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The  remaining  dose-rate  contribution  calculated  with  this  program  is  that  due  to 
the  inscatter  of  gamma  rays  at  the  corner  lip,  as  shown  in  Figure  5.  The  equation 
developed  in  Reference  1  for  this  dose-rate  contribution  is 


where  ri ,  t2»  Sq/  9#  and  6^  are  as  shown  in  Figure  5  and  are  measured  at  the  tip 
of  the  corner  lip;  H  is  the  height  of  the  duct;  RL  is  the  mean  free  path  of  the  gamma 
rays  in  concrete  and  is  given  in  Figure  4;  and  ZN  is  the  number  of  electrons  per  unit 
volume  and  was  calculated  to  be 

ZN  =  1.15  X  10^^  electrons/in. ^ 

3 

for  145-lb/ft  concrete. 

The  dose  rates  calculated  with  this  program  are  about  one-third  of  measured  dose 
rates  for  small  ducts  and  about  one-half  for  large  ducts.  This  would  be  expected  from 
experimental  evidence  showing  the  importance  of  multiple  scattering  in  a  duct. 7 


source 


Figure  5.  Geometry  showing  inscatter  of  gamma  rays 
from  corner  lip. 
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PROGRAM  2.  CALCULATION  OF  PRIMARY  SCATTERING  AREAS  BY  INCREMENTS 


In  the  LeDoux-Chilton  program,  r^ ,  r2,  0q,  0,  and  cp  were  considered  constant 
over  each  primary  scattering  area.  This  would  not  be  a  bad  approximation  for  ducts 
with  very  long  legs,  but  could  give  erroneous  answers  for  ducts  with  legs  of  medium  or 
short  lengths.  Therefore  a  program  was  written  in  which  each  primary  scattering  area 
is  divided  into  small  increments  for  its  calculation.  The  corner  is  divided  into  the 
five  areas  shown  in  Figure  6.  The  ceiling  (or  floor)  is  divided  into  three  areas  for 
ease  of  calculation.  Each  area  shown  in  Figure  6  is  divided  into  several  small  Increments. 
The  dose  rates  for  all  the  increments  are  calculated  using  Equation  1  and  summed  to 
give  the  dose  for  the  primary  scattering  area  under  consideration.  For  each  calculation, 
r] ,  r2,  9q,  0,  and  <p  are  measured  at  the  center  of  the  incremental  area. 

The  dose  rate  from  a  primary  scattering  area  when  calculated  using  small 
incremental  areas  was  generally  50  percent  higher  than  when  calculated  as  one  area, 
and  was  sometimes  as  much  as  100  percent  higher. 


Figure  6.  Division  of  primary  scattering  areas  for  calculation  by 
incremental  areas. 
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For  the  above  results,  eoch  of  the  areas  of  Figure  6  was  divided  into  16  increments. 
It  is  felt  that  further  division  is  undesirable,  because  it  would  greatly  Increase  the 
running  time  on  the  computer.  Also,  finer  divisions  would  probably  give  no  Increase 
in  accuracy  because  of  the  tacit  assumption  that  all  gamma  rays  are  backscattered  at 
the  interface.  Actually,  a  backscattered  gamma  ray  may  emerge  at  distances  of  a 
mean  free  path  or  more  from  the  point  where  it  penetrated  the  wall  of  the  duct. 


PROGRAM  3.  MULTIPLE  SCATTERING  CALCULATIONS 

The  dose-rate  contribution  from  gamma  rays  which  have  been  multiply  scattered, 
that  is  which  have  been  scattered  from  more  than  one  scattering  surface  before  reaching 
the  detector,  has  been  shown  experimentally  to  be  large.  Therefore  it  is  necessary  to 
make  some  attempt  to  calculate  this  dose-rate  contribution. 

A  program  was  written  which  calculates  the  dose  rate  from  gamma  rays  that  have 
been  scattered  from  two  scattering  surfaces.  The  duct  is  divided  into  nine  scattering 
surfaces,  as  shown  in  Figure  7.  All  possible  combinations  are  taken  for  which  a  gamma 
ray  from  the  source  can  scatter  from  one  surface  to  another  surface  and  then  to  the 


Figure  7.  Scattering  areas  used  in  multiple  scattering  program. 
Areas  2,  5,  and  8  are  either  ceiling  or  floor. 
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The  dose  rate  at  the  center  of  the  second  area  due  to  gamma  rays  from  the  source 
impinging  on  the  first  area  is  calculated.  Then  the  dose  rate  at  the  detector  due  to  the 
second  surface  and  the  previously  calculated  initial  dose  rate  is  computed.  The  energy 
of  the  gamma  rays  impinging  on  the  second  area  is  assumed  to  be  the  energy  of  a  gamma 
ray  having  one  Compton  scatter  at  the  center  of  the  first  area  and  going  to  the  center 
ef  the  second  area. 

The  equations  used  for  these  calculations  are  then: 


°  2  2  2 


■■l  ^2  '■3 


26 


C(EQ)K(e^^)10‘''  +  C(EQ) 


cos  6 


1  + 


01 


COS0, 


26 


a«  = 


C(E^)K(0^2)’O 


cos  0, 


1  + 


02 


cos  0. 


where 


A. 

J 


E 

s 


the  area  of  the  first  scattering  surface 
the  area  of  the  second  scattering  surface 

the  energy  of  a  gamma  ray  undergoing  one  Compton  scatter  in  going 
from  the  source  to  the  center  of  A;  to  the  center  of  Aj 


and  r^ ,  x^,  rg,  ,  0q2<  0^ ,  00/  <P]r  and  <^2  are  as  shown  in  Figure  8  for  the 
case  of  i  =  3  and  J  =  5.  C  and  C' are  calculated  from  Equations  3  and  4. 


Three  possible  combinations  of  i's  and  J's  were  found  to  be  too  complicated  to 
be  put  into  this  program.  These  are;  i  =  1  and  j  =  8,  i  =  2  and  j  =  7,  and  i  =  2 
and  j  =  8.  Those  combinations,  believed  to  be  small,  were  therefore  neglected. 
This  lack  is  offset  by  the  assumptions  that  surface  8  is  in  view  of  all  of  surface  4, 
surface  6  is  in  view  of  all  of  surface  2,  and  surface  5  is  in  view  of  all  of  surface  2. 
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Figure  8.  Geometry  for  multiple  scattering  for  i  =  3  and  j  »  5. 

Special  means  were  needed  to  calculate  the  combination  I  =  1  and  j  =  7. 
From  Figure  9; 


Wj 

W  L _ - 

=  (S  -  X,)(Lj  -  -J  -  Xj)  *  W,Wj  In  (-4-^) 

-  >«2(S  -  ^-S) 
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where  A1A7  is  the  effective  product  of  those  parts  of  A-]  and  A7  that  are  in  view  of 
each  other.  For  this  case  r^ ,  ^2,  r3,  0q^,  0p2,  0],  02/  and  0^2  are  not 
measured  at  the  center  of  Ai  and  A7,  but  as  shown  in  Figure  9. 

Dose  rates  calculated  with  the  multiple  scattering  program  generally  represent 
a  sizable  portion  of  the  total  dose  rate.  For  6-foot-square  ducts  calculated,  multiple 
scattering  doses  are  about  one-third  of  the  total  calculated  dose  rate. 


Figure  9.  Geometry  for  calculating  the  1-to-7  contribution. 


PROGRAM  4.  MULTIPLE  CORNER  INSCATTERING  CALCULATIONS 

The  need  for  a  more  inclusive  calculation  of  the  corner  lip  effect  was  shown 
when  the  sum  of  the  previous  calculations  (the  corner  lip  penetration  and  inscatter 
from  the  LeDoux-Chilton  program,  the  primary  scattering  surfaces  from  incremental 
areas,  and  the  multiple  scattering)  gave  fairly  good  results  for  large  ducts,  but  very 
low  results  for  small  ducts.  It  was  believed  that  the  LeDoux-Chilton  method  of 
treating  the  corner  lip  was  correct,  but  that  an  important  contribution  had  not  been 
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considered.  The  importance  has  already  been  shown  for  those  dose  contributions 
involving  two  processes;  i.  e.,  two  backscatters,  and  a  backscatter  and  a  penetration 
of  the  corner  lip  (as  for  A5,  A^,  k-j,  and  Ag  in  the  LeDoux-Chilton  analysis). 
Therefore,  it  was  believed  necessary  to  calculate  those  interactions  involving  back¬ 
scatter  from  a  surface  and  an  inscatter  from  the  corner  lip. 

The  dose-rate  contribution  of  the  four  cases  involving  these  processes  are 
calculated  in  this  program.  These  four  cases  are  shown  in  Figure  10.  The  first 
two  cases  involve  a  backscatter  from  a  surface  in  the  first  leg  to  the  corner  lip 
and  an  inscatter  from  the  corner  lip  to  the  detector.  The  other  two  cases  involve 
an  inscatter  from  the  corner  lip  to  a  surface  in  the  second  leg  and  a  backscatter 
from  that  surface  to  the  detector. 

The  equations  used  for  these  calculations  are,  for  Cases  1  and  2: 


D 


D^ZNK{E,,  ©52^  cos 0^  cosejCRL^)^ 

■  2 

^3 


where 


DpaAj  COS0Q 
2'-2 
’’l  '‘2 


and  for  Cases  3  and  4: 


D  = 


O.aA.  005  0, 

J.  A 


r 


3 


where 


= 


DqZNK(Eq,  e^,)Hcos0^  cos02(RLq)" 

2  “2 
‘‘l  ’’2 
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In  fhe  above  equations: 


a  =  the  albedo,  calculated  by  Equations  2,  3,  and  4 
A.  =  the  area  of  the  scattering  surface  in  the  first  leg 
Aj  =  the  area  of  the  scattering  surface  in  the  second  leg 
D  =  the  dose  rate  at  the  detector 
Dq  =  the  dose  rate  at  one  unit  length  from  the  source 
Eq  =  the  initiol  energy  of  gamma  rays  from  the  source 
E^  =  the  energy  of  a  gamma  ray  scattered  through  6^^ 

K(E,  0  )  =  the  Klein-Nishina  scattering  coefficient 
s 

RLq  =  the  mean  free  path  for  gamma  rays  of  energy  Eq 

RL^  =  the  mean  free  path  for  gamma  rays  of  energy  E^ 

ZN  =  the  number  of  electrons  per  unit  volume  in  the  corner  lip 

and  ri ,  ^21  rg,  9q,  0,  0i ,  So,  65I ,  tp  are  as  shown  in  Figure  10  and 

ore  measured  at  the  center  of  the  scattering  area  under  consideration  and  at  the  tip 
of  the  corner  lip. 

These  calculations  give  large  dose-rate  contributions,  as  will  be  seen  in  the 
following  analysis  of  results.  In  general  they  account  for  about  50  percent  of  the 
total  dose  rote  in  small  (1 1 -Inch-square)  ducts  and  about  15  percent  in  large  (6-foot- 
square)  ducts.  The  contributions  from  Cases  3  and  4  are  5  to  10  times  the  contributions 
from  Cases  1  and  2. 


COMPARISON  OF  CALCULATIONS  AND  EXPERIMENTAL  DATA 

The  results  of  oil  the  calculations  performed  to  date  are  given  in  Table  I.  As 
can  be  seen,  the  agreement  between  calculated  values  and  measured  values  is  normally 
better  than  ±30  percent.  The  exceptions  to  this  agreement  are  the  12-inch  duct  with 
the  Cs-137  source,  ducts  with  very  short  first  legs  (L]/W  =  1.33),  and  the  values  for 
Au-198  sources.  This  is  considered  good  agreement,  as  error  in  experimental  measure¬ 
ments  and  source  strengths  would  probably  introduce  errors  of  10  to  15  percent.  This 
can  be  seen  by  comparing  the  Co-60  data  of  Green, ^  Eisenhauer,  ®  and  Terrell  ^  for 
Li/W  ~  3.5. 
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Figure  10.  Geomelry  for  multiple  corner  inscattering  with  scattering  surface  in  the 
first  leg  (Cases  1  and  2)  and  in  the  second  leg  (Cases  3  and  4). 

Areas  2  and  5  shown  in  Case  1  are  either  ceiling  or  floor. 
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The  larger  error  for  the  12-inch  duct  with  the  Cs-137  source  is  evidently  due 
to  an  overestimation  of  the  corner  lip  effect.  The  use  of  the  energy-absorption 
coefficient  to  obtain  the  mean  free  path  is  an  overestimation  that  might  become 
important  in  small  ducts  at  lower  incident  gamma-ray  energies.  Also,  since  the 
energy-absorption  coefficient  is  used,  there  is  some  double  counting  of  gamma  rays 
from  A5 ,  A^,  Ay,  and  A0  by  including  these  areas  in  the  geometry  for  the  multiple 
corner  lip  inscattering  program. 

The  fact  that  the  calculations  are  more  in  error  for  very  short  first  legs  can  be 
attributed  to  the  assumptions  going  into  the  calculations;  mainly  those  simplifying 
the  areas  in  the  multiple  scattering  program,  and  the  assumption  that  the  gamma  rays 
impinging  on  the  corner  lip  are  parallel. 

The  discrepancy  for  Au-198  is  more  difficult  to  explain  on  the  basis  of  the 
calculations.  The  fact  that  there  is  something  peculiar  about  the  Au-198  measurements 
can  be  seen  from  Figure  1 1. 


0.2  0.3  04  0.5  04  0.8  1  3 

Gamma-Ray  Energy  (mev) 


Figure  1 1.  Variation  of  dose  attenuation  with  energy  for  a  given  duct. 
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This  figure  is  essentially  the  graph  of  experimentally  determined  attenuation 
factors  in  a  6-foot-square  duct,  with  Li/W  =  2  and  L2/W  =  3.  17,  versus  gamma-ray 
energy.  The  points  for  Na-24  and  Au-198  were  obtained  by  interpolation  of  Terrell's 
data.  The  point  for  Au-198  is  much  lower  than  one  would  expect  from  the  trend  of 
the  curve.  The  peculiarity  of  this  point  could  be  due  to  natural  phenomena,  or  to 
error  in  the  measurements.  More  experimental  measurements  of  Au-198  should  be 
made. 


Tables  II  through  V  indicate  the  general  magnitudes  of  the  many  contributing 
sources  of  scattered  radiation  in  o  duct.  The  contributions  from  all  of  the  scattering 
sources  considered  in  these  calculations  have  been  listed  in  Tables  II  and  III  for  Co-60 
and  Tables  IV  and  V  for  Cs-137,  both  for  a  small  and  a  large  duct. 

The  dose-rate  contributions  calculated  for  multiple  corner  lip  scattering  are 
very  large,  especially  forsmoll  ducts.  That  the  contribution  from  the  corner  lip 
is  large  in  actuality  can  be  seen  from  Figure  12. 


Wid»h,  W  (in.) 

Figure  12.  Variation  of  corner  lip  effect  with  width  of  duct. 
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In  this  figure  the  attenuation  coefficient,  DL^^/Dq,  for  L]/W  =  L2/W  =  2  is 
plotted  versus  W.  Also  plotted  in  this  figure  is  the  calculated  (Dp^  +  Li^/Dq, 
w/here  Dp^  is  the  dose  from  the  primary  scattering  areas  and  is  the  dose  from 
multiple  surface  scattering.  Since  (Dp^  +  D/y/\s)/Do  scales  as  1/L^^,  the  calculated 
(Dpy^  +  D^5)L]2/Dq  is  constant  on  the  graph,  and  since  the  only  assumption  is  that 
scattered  gamma  rays  enter  and  exit  at  the  same  point,  one  \would  expect  that  the 
physical  values  for  (Dp^  +  /Dq  would  also  plot  as  a  constant  if  they  could 

be  found.  Since  the  corner  lip  effect  decreases  with  W  (about  as  l/W),  the  experi¬ 
mental  DL]2/Dq  does  appear  to  asymptotically  approach  a  constant,  which  is  drawn 
on  the  figure.  This  constant  is  therefore  labeled  the  experimental  (Dp  a  +  D|^5)Li2/Do/ 
and  any  dose  over  this  constant  is  due  to  the  corner  lip  effect.  As  can  ee  seen,  tne 
corner  lip  effects  are  by  far  the  major  contributors  for  small  ducts. 

The  multiple  surface  scattering  also  gives  a  large  contribution.  As  can  be  seen 
in  the  breakdown  of  multiple  scattering  contributors,  this  is  due  to  24  small  contributors 
adding  up  to  a  sizable  value,  rather  than  to  a  few  large  contributors. 


CONCLUSION 

In  summary,  these  programs  appear  to  compute  the  dose  rates  in  ducts  through 
concrete  with  sufficient  accuracy  to  assist  in  entranceway  design.  The  calculations 
also  point  the  way  to  fruitful  experimentation,  such  as  resolving  the  Au-198  anomaly 
and  determining  the  actual  corner  lip  effect  in  small  ducts. 
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Table  I.  Comparison  of  Calculated  and  Measured  Dose  Rates 


Difference  =  (measured 'coicuioted/meowrtd)  100. 


2/  The  volues  for  Na'24  are  the  sum  of  the  values  obtoined  using  initio!  gommo-roy  energies  of  1. 37  mev 
and  2. 75  mev. 


19 


6  5 


Scattering  %  of 
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Scattering 
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1 - 

%  of  Multiple 
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2-3 
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3-5 

55.0 
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99.9 

inscatter 

Multiple  19.0 

surface 

scatter 

Multiple  45.7 
corner  lip 
scatter 


100.1 
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Table  III.  Comparison  of  Scattering  Source  Contributions  in  o  6-Foot  Duct: 
L^/W  =  2,  Lj/W  =  2.5;  Co-60  Source 

4  1  1 


Scattering 

Source 

%  of 
Total 

A] 

A2 

12.8 

11.2 

A3  + A4 

23.5 

A5 

0.8 

A^ 

1.0 

4 

0.3 

As 

1.7 

Single 
corner  lip 
Inscatter 

3,8 

Multiple 

surface 

scatter 

27.8 

Multiple 
corner  lip 
inscatter 

16.9 

Total 

99.8 

Breakdown  of  Multiple 
Scatter  Contribution 


Scattering 

Combination 


%  of  Multiple 
Scatter 
Contribution 


Breakdown  of  Multiple 
Corner  Lip  Contribution 

Scattering 

Combination 

%  of  Multiple 
Corner 
Contribution 

1-3 

9.0 

2-3 

8.7 

3-4 

32.5 

3-5 

49.8 

Total 

100.0 

100.2 
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Table  V.  Comparison  of  Scatfering  Source  Contributions  in  a  6-Foot  Duct: 
L^/W  =  2,  L2/W  =  2.5;  Cs-137  Source 


Breakdown  of  Multiple 
Scatter  Contribution 

Scattering 

Combination 

%  of  Multiple 
Scatter 
Contribution 

1-5 

5.7 

1-6 

2.2 

1-7 

5.7 

2-4 

6.6 

2-5 

5.5 

2-6 

4.3 

3-4 

3.2 

3-5 

4.5 

3-6 

2.1 

4-5 

3.6 

4-6 

3.4 

4-7 

4-8 

4-9 

5-4 

5-5 

1.3 

5-6 

3.8 

5-7 

6. 1 

5-8 

6.8 

5-9 

6-4 

6-5' 

6-8 

6. 1 

6-9 

Total 

—■■—Mil 

Scattering 

Source 

%  of 
Total 

12.7 

11.0 

A3  + A4 

22.0 

0.9 

4 

0.9 

A7 

0.3 

^8 

1.4 

Single 
corner  lip 
inscatter 

3.9 

Multiple 

surface 

scatter 

28.5 

Multiple 
corner  lip 
inscatter 

18.4 

Total 

100.0 

Breakdown  of  Multiple 
Corner  Lip  Contribution 

Scattering 

Combination 

%  of  Multiple 
Corner 
Contribution 

1-3 

8.9 

2-3 

8.3 

3-4 

31. 1 

3-5 

51.7 

Total 

100.0 

23 


REFERENCES 


1.  NCEL  Technical  Note  N-383,  Attenuation  of  Gamma  Radiation  Through 
Two-Legged  Rectangular  Ducts  —  An  Analytical  Approach,  by  J.  C.  LeDoux  and 
A.  B.  Chilton.  Port  Hueneme,  California,  20  January  1961. 

2.  Technical  Operations,  Inc.  Report  TO-B  61-39,  Monte  Carlo  Calculations  on 

the  Reflection  and  Transmission  of  Scattered  Gamma  Radiation  (revised),  by  D.  J.  Raso. 
Burlington,  Massachusetts,  May  1962. 

3.  NCEL  Technical  Report  R-228,  A  Semiempirical  Formula  for  Differential  Dose 
Albedo  of  Gamma  Rays  on  Concrete,  by  A.  B.  Chilton  and  C.  M.  Huddleston. 

Port  Hueneme,  California,  16  November  1962. 

4.  NCEL  Technical  Note  N-469,  Some  Applications  of  a  Semiempirical  Formula 
for  Differential  Dose  Albedo  for  Gamma  Rays  on  Concrete,  by  W.  C,  Ingold. 

Port  Hueneme,  California,  30  November  1962. 

5.  Rockwell,  T.,  Ill,  editor.  Reactor  Shielding  Design  Manual.  McGraw-Hill, 

New  York,  1956. 

6.  NCEL  Technical  Report  R-195,  Attenuation  of  Gamma  Radiation  in  a  Two-Legged 
ll-lnch  Rectangular  Duct,  by  D.  W.  Green.  Port  Hueneme,  California,  2  May  1962. 

7.  NCEL  Technical  Note  N-443,  Gamma  Dose  Rates  and  Energy  Spectra  in  a  3-Foot 
Square  Duct,  by  J.  M.  Chapman.  Port  Hueneme,  California,  29  June  1962. 

8.  National  Bureau  of  Standards.  Technical  Note  74,  Scattering  of  Cobalt-60 
Gamma  Radiation  in  Air  Ducts,  by  C.  Eisenhauer.  Washington,  D.  C. ,  1960. 

9.  Armour  Research  Foundation.  Report  ARF  1 158-A02-7,  Radiation  Streaming  in 
Shelter  Entranceways,  by  C.  W,  Terrell,  et  al.  Chicago,  1962. 

10.  Armour  Research  Foundation.  Report  ARF  1 158-12,  Radiation  Streaming  in 
Shelter  Entranceways,  by  C.  W.  Terrell,  et  al.  Chicago,  1960. 

11.  Armour  Research  Foundation.  Report  ARF  1158-A01-5,  Radiation  Streaming  in 
Shelter  Entranceways,  by  C.  W.  Terrell,  etal.  Chicago,  1961. 


24 


DISTRIBUTION  LIST 


SNDL 

No .  of 

Total 

Code 

Activities 

Copies 

1 

10 

Chief,  Bureau  of  Yards  ond  Docks  (Code  42) 

23A 

1 

1 

Naval  Forces  Commanders  (Taiwan  only) 

39  B 

2 

4 

Construction  Bottolions 

39  D 

10 

10 

Mob!  le  Construction  Battalions 

39  E 

3 

3 

Amphibious  Construction  Battolions 

39  F 

1 

2 

Construction  Battalion  Base  Units 

A2A 

1 

1 

Chief  of  Navol  Reseorch  •  Only 

A3 

2 

2 

Chief  of  Novol  Operation  (OP-07,  OP-04) 

AS 

S 

5 

Bureaus 

B3 

2 

2 

Colleges 

E4 

1 

2 

Loboratory  ONR  (Washington,  D.  C.  only) 

E5 

1 

1 

Research  Office  ONR  (Pasadena  only) 

E16 

1 

1 

Training  Device  Center 

F9 

7 

7 

Station  ‘  CNO  (Boston;  Key  West;  Son  Juan;  Long  Beach;  San  Diego; 

T reasure  Island;  and  Rodman,  C .  Z  .  only) 

FI7 

6 

6 

Communication  Station  (San  Juan;  Son  Francisco;  Pearl  Harbor;  Adok, 
Alaska;  and  Guam  only) 

F41 

1 

1 

Security  Stotlon 

F42 

1 

1 

Rodio  Station  (Oso  and  Cheltanham  only) 

F48 

1 

1 

Security  Group  Activities  (Winter  Harbor  only) 

H3 

8 

8 

Hospitoi  (Chelsea;  St.  Albans,  Portsmouth,  Va;  Beaufort;  Great  Lokes; 
Son  Diego;  Oakland;  ond  Camp  Pendleton  only) 

H6 

1 

1 

Medical  Center 

J1 

2 

2 

Adminiatrotion  Commend  ond  Unit  -  BuPers  (Great  Lakes  ond 

San  Diego  only) 

J3 

1 

1 

U.  $.  Fleet  Anti-Air  Worfore  Troining  Center  (Virginia  Beoch  only) 

J4 

2 

2 

Amphibious  Boses 

J19 

1 

1 

Receiving  Stotion  (Brooklyn  only) 

J34 

I 

1 

Stotlon  •  BuPers  (Washington,  D .  C  .  only) 

J46 

1 

1 

Personnel  Center 

J48 

1 

1 

Construction  Troining  Unit 

J60 

1 

1 

School  Academy 

J65 

1 

1 

School  CEC  Officers 

J84 

1 

I 

School  Postgroduote 

J90 

1 

1 

School  Supply  Corps 

25 


DISTRIBUTION  LIST  (Cont’d) 


SNDL 

No .  of 

Total 

Code 

Activities 

Copies 

J95 

1 

1 

School  War  College 

J99 

1 

1 

Communication  Training  Center 

LI 

11 

11 

Shipyards 

L7 

4 

4 

Laboratory  -  BuShips  (New  London;  Panomo  City;  Corderock;  and 
Annapolis  only) 

L26 

5 

5 

Naval  Facilities  -  BuShips  (Antiguo;  Turks  Island;  Barbados; 

Son  Salvador;  ond  Eleuthera  only) 

L30 

1 

1 

Submarine  Base  (Groton,  Conn,  only) 

L32 

2 

2 

Naval  Support  Activities  (London  and  Naples  only) 

L42 

2 

2 

Fleet  Activities  >  BuShips 

M27 

4 

4 

Supply  Center 

M28 

6 

6 

Supply  Depot  (except  Guantanamo  Bay;  Subic  Bay;  end  Yokosuka) 

M61 

2 

2 

Aviotion  Supply  Office 

N1 

6 

18 

BuDocks  Director,  Overseas  Division 

N2 

9 

27 

Public  Works  Offices 

NS 

3 

9 

Constructior>  Battalion  Center 

N6 

5 

S 

Construction  Officer>in*Charge 

N7 

1 

1 

Cert  struct  ion  R  esi  dent ‘Off  icer'in^Charge 

N9 

6 

12 

Public  Works  Center 

N14 

1 

1 

Housing  Activity 

R9 

2 

2 

Recruit  Depots 

RIO 

2 

2 

Supply  Installotions  (Albany  and  Borstow  only) 

R20 

1 

1 

Marino  Corps  Schools  (Quantico) 

R64 

3 

3 

Marine  Corps  Base 

R66 

1 

1 

Marine  Corps  Camp  Detochment  (Tengan  only) 

WlAl 

6 

6 

Air  Station 

W1A2 

35 

35 

Air  Station 

WIB 

8 

8 

Air  Station  Auxiliary 

W1C 

4 

4 

Air  Facility  (Phoenix;  Monterey;  Oppoma;  Naha;  and  Naples  only) 

W1E 

6 

6 

Marine  Corps  Air  Station  (except  Ouontico) 

WIH 

9 

9 

Station  •  BgWeps  (except  Rota) 

1 

1 

Deputy  Chief  of  Staff,  Reseorch  and  Development,  Headquorrers, 

U  .  S  .  Mo  fine  Corps,  Wash  ing  ton,  D .  C  • 

1 

1 

President,  Morine  Corps  Equipment  Board,  Marine  Corps  School, 

Quonti  CO,  Va  . 

1 

2 

Library  of  Congress,  Washington,  D.  C. 

1 

1 

U.S,  Naval  Research  Laboratory,  Chemistry  Division,  Washington,  D.C. 

1 

1 

Deputy  Chief  of  Staff,  Research  and  Development  Headquarters, 

U.S«  Marine  Corps,  Washington,  D.C. 

26 


No.  of 
Acti  viti  es 


DISTRIBUTION  LIST  (Cont’d) 


1 

1 

1 

1 


Total 
Copi  es 

1 


1 

1 

1 

3 

1 

1 

1 

2 

20 

2 

2 

2 

1 

1 

2 

I 

1 

1 

I 

I 

4 

1 

1 

1 

1 


Chief  of  Staff,  U.  S.  Army,  Chief  of  Research  and  Development,  Department  of 
the  Army,  Washington,  D.  C« 

Office  of  the  Chief  of  Engineers,  Assistant  Chief  of  Engineering  for  Civil  Works, 
Department  of  the  Army,  Washington,  D.  C. 

Chief  of  Engineers,  Department  of  the  Army,  Washington,  D  ,  C .,  Attn:  Engineering 
Research  and  Development  Division 

Chief  of  Engineers,  Department  of  the  Army,  Washington,  D.  C.,  Attn:  ENGCW-OE 

Director,  U.S.Army  Engineer  Research  and  Development  Laboratories, 

Fort  Belvoir,  Va.,  Attn:  Information  Resources  Branch 

ASD  (ASNXRR),  Wright-Potterson  Air  Force  Base,  Ohio 

Headquarters,  U.S.Air  Force,  Directorate  of  Civi  I  Engineering,  Washington,  D .  C 
Attn:  AFOCE-ES 

Commanding  Officer,  U.S,  Naval  Construction  Battalion  Center,  Port  Hueneme, 
Calif.,  Attn;  Materiel  Dept.,  Code  140 

Deputy  Chief  of  Staff,  Development,  Director  of  Reseorch  end  Development, 
Department  of  the  Air  Force,  Washington,  D.  C. 

Director,  National  Bureau  of  Standards,  Department  of  Commerce,  Connecticut 
Avenue,  Washington,  D.C. 


Office  of  the  Director,  U,  S.  Coast  and  Geodetic  Survey,  Washington,  D.  C* 

Defense  Documentation  Center,  Building  5,  Comeron  Station,  Alexondrio,  Va  , 

Director  of  Defense  Research  ond  Engineering,  Department  of  Defense, 
Washington,  D.  C. 

Diteetor,  Division  of  Plons  and  Policies,  Headquorters,  U.S.  Marine  Corps, 
Washington,  D.C. 

Director,  Bureau  of  Reelomatien,  Washington,  D.  C. 


Facilities  Officer,  Code  108,  Office  of  Naval  Research,  Washington,  D.  C. 

Federal  Aviation  Agency,  Office  of  Monogement  Services,  Administrative  Services 
Division,  Washington,  D.  C.,  Attn:  Library  Bronch 


Norfolh^'^Va.^'’'''’'  Amphibious  Bose,  Little  Creek, 


PornTMufl":  documentation  Section,  P.O.  Bo*  10, 

Air  Force  Cambridge  Research  Center,  Hanscom  Field,  Bedford,  Mass. 

Directorate  of  Research,  Air  Force  Special  Weopons  Center,  Kirtland  Air  Force  Bose,  N.M. 
Sondia  Corporotion,  Box  5800,  Albuquerque,  N.M,  Attn:  Classified  Document  Division 
Director,  Engineering  Research  Institute,  University  of  Michigan,  Ann  Arbor,  Mich. 

U.S.  Army  Material  Command,  Washington,  D.C. 

Dr.  H.  E.  Stanton,  Physics  Division,  Argonne ’National  Laboratory,  Argonne,  111. 

Mr,  John  H.  Hubbell,  Notional  Bureau  of  Stondards,  Washington,  D.C. 

Professor  D.  W.  Green,  1668  Maple  Avenue,  Galesburg,  III. 

Rivers  and  Harbor  Library,  Princeton  University,  Princeton,  N.  J. 


27 


§s 


o 

CN 
CO 

^  >  I 

*rt  to 

C  c9 

Sol 

.2  i-o 

-o  o  lY 

n  JC  • 

Su  >■ 


o 

O 

i 


UJ  O 

<  2 

O'  S 

iA  < 
O  UJ 

a  X 


cn 

»/> 

< 

U 

z 

3 


o 

s 

u  - 


$ 

-o 


f  8. 


O)  0 

4)  O 

C  ^ 

■§>3  < 

~  “u 

^  8.S 


n  hi 


O  § 

2  1.3 

U1  0  4- 

3  J=  U 
U  O 

u  2 

3 

Q  3 

uu  X 

O  , 
a  £ 

UJ  lu 

-I'  y  -. 
O  Z  a. 
>  O 
^  U 


2  s£ 


■S  o 


6  •?  S 


4>  G  . 


lA  -5 

.2  a. 


6  £ 
U  = 


S  o 

ij 

01  g- 

(U  X 

Jl  4) 


0) 

s 

a 

I 

t: 


CM 


o  Q.  o 

-o  OU. 

n  -C  I 

i2  u  >- 


o 


i 


O 

Q 


o , 

a 

u.  ' 

O 

z 

o 


u 

z 


• - D 

M  2 


■§>5< 


t  ^ 

^  8.5 

—  u 
o  £K  q: 

>  _  UJ 

■7  °  3 

cl  o  ~ 


.c  < 

'J  o 

3  w  U 


O  § 

z  I.  2 

uo  a  4- 
3  jC  u 
° 

u  2  « 
3 

Q  -; 

Q  X 

UJ  J3 

o  . 
ot  S 

UJ  UJ  = 

■f  U  ••; 
0  2  o- 
>  O  1^ 
^  u 


S. 


4) 

p  — 

xi 

0  u 


3i  p  6 

-S 

g  dj  2 
.9  ‘5  2; 
0  4>  2^ 

-ii  i  -u 

o  -  c 
y  ^  O 

i.  0  •» 

<u  jii 

*“  o  ^ 
a*  o)  ^ 

E  -Q  T3 

a  o  ^3 
U  4> 

4>  G> 

S  D3 

O  ns 

“D  — 


2.S 

4)  i  ^  On 

i 

'  Ic^ 

e  /§  0  ^ 

^  .9  iz 

4)  "O  O  'Y 

iS  tS  j5  >■ 


o 

e 

o 

JO 

o 


a  0 

UJ  3! 

U1  < 
O  UJ 

a  X 

U.  •" 

O  O 


to 

in 

< 

3 

(J 

Z 

3 


ii 
£  2 
^  c 


4) 

X  P 


23 


0 

E  CO 
Q.  «0 

-8  t 
u  o 


■D 

«J 

u 

8. 

E 

O 

S 

o 


.Sex 

n  — 


LU  «N 


3!^  2 


> 

u 


u 


Q 


8.5° 


1.U 

;  a;  O 

1  3 


Z  .  _ 

•  =  %  A 
u^  .c  S  O 


—4  u 

ID  H“ 


a  o  § 


UJ  3 

u  - 

z  d 

O  X 
(J  CN 


x:  E 
8  SI 

u  O 

i  o 
8.-8 
E  :S 
6  ° 
s 

O 

T) 


5  C 

"I 

«>  X 

J:  4) 


s 

4> 

s 

Q 

I 


E  ^  O 

□  <  >Ci 


is! 

y4-|r- 

j:  ,2  Q-  o 

4)  ■G  2  ^ 

tS  u  > 


UJ  •-  I 
—  t;  ft 


^ 

0 

p 


Su  o 

t;  a 


UJ 


It 


u 

z 

3 


-  TJ 

^  $  u 

£  D  8. 

u  u  g 
C  O 

o  P  O 


is 

UJ  ^ 


O  O  g 
zZ  I 

n  tiO  o 

-  =>  7^ 


■D 


P 


^  I 

ri  — 


u  o 


t 


u 


-  <  o 

a.^  UJ 

^  O 
t;  ^ 

a  t:  UJ 
a  3  Zj 


5  cs  -li  -5  ^ 

^  2  o  *“ 

♦  £  T3 


w  - 

J>.  ^ 


ui‘  A  2  o 

3  I-  U  I— 


U  -. 
Z  a 
O  h- 
u 


E  £ 
o 


U 


il 


:perimentol  data.  t  {  experimental  data. 


